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d 1. Introduction 



ABSTRACT 



Context. Light elements and nitrogen surface abundances together can constrain the mixing efficiencies in massive stars on the main 

sequence, because moderate mixing in the surface layers leads to depletion of light elements but only later to enrichment in nitrogen. 

Aims. We want to test the rotational mixing prescriptions included in the Geneva stellar evolution code (GENEC) by following the 

evolution of surface abundances of light isotopes in massive stars. 

Methods. The GENEC is a ID code containing sophisticated prescriptions for rotational mixing. We implemented an extended 

reaction network into this code including the light elements Li, Be, and B, which allowed us to perform calculations testing the 

rotation-induced mixing. 

Results. We followed 9, 12, and 15 M© models with rotation from the zero age main sequence up to the end of He burning. The 

calculations show the expected behaviour with faster depletion of light isotopes for faster rotating stars and more massive stars. 

Conclusions. We find that the mixing prescriptions used in the present rotating models for massive single stars can account for most 

of the observations; however, the uncertainties are quite large, making it hard to draw a firm conclusion on the mixing scenario. 

Key words. Stars: abundances - Stars: massive - Stars: rotation 
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Rotation is beside the stellar mass and the initial chemical com- 
position a key parameter in the evolution of single stars. It af- 
fects the physical and chemical structures of the stars and there- 
fore quantities such as lifetime, luminosity, effective tempera- 
ture Teff etc. Recent models including rotation reproduce a wide 
range of observations better than those without as for example, 
the rutrogen surface enrichment ( Heger & Langer 2000 ; Meynet 
|& Maeder] |2000[), th e Wolf-Ray et to O-type star numb er ratio 
( [Meynet & Maederj |2003 , 2005; [Vazquez et aE] [2007] l, the ra- 
tio of blue to red supergiants in the Small Magellanic Cloud 
( [Maeder & Meynet| [200 l| l, or the variation with the metallic- 
ity of the number ratio of type Ibc to type II supernovae ( Georgy 
[et al. 2009 1. Still, the treatment of transport of angular momen- 
tum and chemical species is thought to be one of the main un- 
certainties in stellar evolution models. New observational data of 
late O- and B-type stars from the VLT-FLAMES survey ([Evans 



Therefore modest mixing due to rotation leads to a depletion of 
light elements at the surface without considerable nitrogen en- 
richment. This effect can not be explained by mass transfer in 
a binary system, since there the accreted material is depleted in 
boron and enriched in nitrogen (,Fhegner et al. 1996 1. 



Boron is produced in spallation process of CNO atoms in the 
interstellar medium (ISM) by galactic cosmic rays (GCRs). In 
massive stars boron is only destroyed. It is the only light element 
out of Li, Be and B which is observed at the surface of massive 
stars (OB-type). Despite the difficulties to measure boron sur- 
face abundances an increasing number of boron surface abun- 
dances from O- and early B-type stars became available in the 
last few years ([Mendel eraL|[2006| [Venn et al.[ [2002[ [Proffitt & 



'Qui gley[[2001[[Proffitt eraLJ(l999 V The comparison in 'Mendel 
et al.[ ([2006[) of observational data with the models of Heger 



[etaL][2005i [20061 [Hunter eTaL|[2007i [Trundle et al.||2007| l lead 
to an intense discussion about how well the models with rotation 



can explain the observed nitrogen surface abundances ( [Hunter 
[eTal 2008 2009a. M aeder et al.| [2008l l and about whether or 
not binaries are needed to explain some groups of the observed 
O- and B-type stars (Langer et al. 2008 1. Light elements and 



in particular boron can constrain the mixing induced by rotation 
and help distinguish between single stars and interacting binaries 
( Brott et aLl p009). Boron is destroyed at relatively low temper- 
atures (* 6 ■ 10" K) where the CNO-cycles are not yet efficient. 



[& Langer ( 2000| l shows a good agreement with the exception 
of three stars (HD 30836, HD 36591, HD205021). The strong 
boron depletion in these young stars raises the question if the 
efficiency of surface mixing due to rotation should be stronger 
or if there are other mixing processes at work. We reexamined 
this question because the Geneva stellar evolution code (here- 
after GENEC) includes the effect of rotation in a different way 
with respect to the codes which were used in previous works 
to examine that question. The most important difference comes 
from the fact that in GENEC the transport of the angular mo- 
mentum is properly accounted for as an advection process and 
not as a diffusion process. 
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In this paper we present correlations of the surface boron 



abundances with nitrogen as in Mendel et al. ( 2006 1, Venn et al. 



(|200 2 ) and also with other interesting quantities such as the 
C/'-^C number ratio or observable physical quantities such as 
the gravity and the surface velocity. In Sect.l2]we give a short de- 
scription of the model ingredients and present the set of simula- 
tions performed. In Sect.|3] we discuss the results from our mod- 
els. In Sect. |4]we compare them to observations and in Sect. l5] 
we summarise the results. 



2. Stellar model description 

2.1. Rotation induced mixing 

The Geneva code (GENEC) used for the calculation of our mod- 
els is described in detail in many previous publications as for 
example in Hirschi et al. ( 2004 1 and more recently Eggenberger' 
|et al. (2007 ). Since the mixing efficiency is tested by the light 
element surface depletion, we will briefly describe here the mix- 
ing prescription implemented in GENEC. The horizontal trans- 
port of matter is assumed to be much faster than the vertical one, 
which leads to almost constant angular velocity on isobars. This 
in turn enables to describe the stellar structure by shellular rota- 
tion, which allows to keep the stellar structure equations in one 
dimension ( Zahn 1992[ Meynet & Maeder 1997| |. The transport 
of angular momentum in the radiative zones is then described by 



^s(''"k=5;^i(^''"''«)*575iK"-* 
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where p is the density, fi the angular velocity of a shell, Dshear 
the diffusion coefficient due to the vertical shear turbulence (see 
Eq. 4]below) and U{r) the quantity intervening in the expression 
of the radial component of the meridional velocity which is ex- 
pressed by u(r, 6) = U(r) ■ f 2(cos 6). Meridional mixing is an 
advective process. In contrast with diffusive processes which al- 
ways smooth gradients, advection can both build up gradients or 
smooth them. It is therefore important to account for this process 
not as a diffusive process but as an advective one. This has been 
properly done in the present work during the main-sequence 
(MS) phase. After the MS, the impacts of meridional currents 
are much less important because the evolutionary timescales be- 
come shorter. During the post MS, the main effect governing the 
evolution of the angular velocity in the radiative zone is simply 
the local conservation of the angular momentum. 

The transport of chemical composition in the convective core 
is treated as instantaneous. The size of the convective core is 
determined by the Schwartzschild criterion to which an over- 
shooting distance dov, which is set by t/ov - ffov min(//p, rcore), is 
added. We adopted ffov = 0.1 for all calculations. 

The change of the mass fraction, Xj, of a nuclide / due to ro- 
tation induced mixing is described by a diffusion equation (note 
that nuclear burning and mixing are treated separately). 
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where D^ix — Z)eff + Oshear with Dsheai accounting for the verti- 
cal mixing introduced by shear turbulence and Deft the diffusion 
coefficient resulting from the interaction of the strong horizontal 
mixing induced by shear turbulence and the meridional currents 
(see Chaboyer & Zahn 1992[l. The expression of Dgff is given by 



£>eff = 
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For the vertical shear diffusion coefficient we use the expression 
of [Talon & Zahn| ( |1997P which accounts for the effects of hori- 
zontal turbulence: 
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where K is the thermal diffusivity, the diffusion coefficient Dh 
describes horizontal turbulent transport. As a standard, we used 



the Dh derived by Maeder ( 2003 1 



2.2. Nuclear reaction network 

The nuclear reaction network used previously in the Geneva 
stellar evolution code included a limited number of isotopes 
and reaction rates. This prevented the investigation of the evo- 
lution of different isotopes as for example the light isotopes 
lithium, beryllium and boron. We therefore implemented into 
GENEC the Basel reaction network developed originally by 
F.-K. Thielemann, which is more flexible in terms of choice of 
nuclei followed and corresponding reaction rates. This reaction 
network was previously used in a wide range of astrophysical 
nucleosynthesis calculations, e.g. Thielemann & Arnett ( 1985|l, 
Frohlich et al. (20061, etc. The equations describing the abun- 
dance changes and the method how these equations are solved 
are described in |Hix & Thie lemann' (I 1999] l. The reaction rates 
are used in their analytical form in the so-called REACLiB-format 
(see |Rauscher & Thielemannl|2000| l. 

The two reactions which determine the burning timescales 
of H- and He-burning, ''*N(;?, y)'^0 and 3a were taken from 



Imbriani et al. (2005) and Fynbo et al. (2005 1 respectively. The 



analytical fits of these two rates were provided by the JINA rea- 
clib website (groups.nscl.msu.edu/jina/reaclib/db). The rates of 
the ip,a) and (p,y) reactions involved in the different CNO- 
cycles were taken from NACRE (j Angulo et al.| |1999| l. Also all 
the (p, a) and (p, y) reactions on the isotopes of lithium, beryl- 
lium and boron which are responsible for their destruction come 
from this source. The involved y6^-decays are experimental rates 
which can also be found in the JINA reaclib database (under la- 
bel "bet+"). 

We included 43 isotopes from hydrogen up to silicon in the 
reaction network. These isotopes are listed in Table fTl ^Be is 
included implicitly, i.e. assumed to decay instantaneously into 
two a-p articles. 

2.3. Model grid 

The most recent abundance determin ations of boron were based 
on the B III 2065.8 A line strength, fvinn et al.| (|2002]l found 



a plateau of the line strength for Teff between 18'000 and 
29'000 K. Such surface temperatures are reached by MS stars 
between about 9 to 15 Mq. Therefore we chose models with ini- 
tial masses of 9, 12 and 15 Mq to investigate the effect of mixing. 
For each mass, models with different rotational velocities were 
calculated. All models were started from the zero age main se- 
quence (ZAMS) and followed up to the end of He-core burning. 
The time averaged equatorial velocities on the MS phase (i^eq) 
are between to 350 km s"' (see Table|2]i. The models with ro- 
tation were started at ZAMS with a flat angular velocity profile. 
Most of our models were started with a solar like composi- 
tion, since the observations of boron in B-type stars have been 
made for stars in the solar neighbourhood. As initial compo- 
sition we chose X=0.72, Y=0.266 and Z=0.014 with the ele- 
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Fig. 1. Diffusion coefficients and angular velocity versus mass for a 12Mq star with i^ini/wcrit - 0.4 at the start of hydrogen burning 
(a) when central hydrogen mass fraction is 0.48 (b) and 5 ■ 10"^ (c). The area below the convective diffusion coefficient Dconv is grey 
shaded to depict the convective core. The diffusion coefficient describing meridional circulation Dejj is plotted with two different 
lines to distinguish between U{r) < (green dotted line) and U{r) > (green dashed line), i.e. the zones transporting angular 
momentum outwards and inwards respectively. 



Table 1. Isotopes considered in the reaction network and their locity Weq, and the surface mass fractions of p, ^He, ^He, '"B, 
initial abundance in mass fractions. "B.'^C, '^C, and '"^N at different times of the evolution. 



Isotope Mass Fraction I Isotope Mass Fraction 



n 

P 

D 

^He 

''He 

^Li 

'Li 

'Be 

8Be 

'Be 

«B 

"B 

"C 

"c 

[4q 

"N 

14n 

'^o 

I6q 



O.OOOE-00 
7.200E-01 
1.397E-05 
4.415E-05 
2.660E-01 
4.004E-12 
5.689E-11 
O.OOOE-00 
O.OOOE-00 
1.692E-10 
O.OOOE-00 
7.786E-10 
3.465E-09 
O.OOOE-00 
2.283E-03 
2.771E-05 
O.OOOE-00 
O.OOOE-00 
6.588E-04 
2.595E-06 
O.OOOE-00 
5.718E-03 



180 
IVp 
18p 
19p 

20Ne 

2'Ne 

^^Ne 

2'Na 

22Na 

23Na 

2*Mg 

^^Mg 

2S]y[g 

25 Al 

2SA1 

2'Al 

"Si 

2«Si 

^'Si 

^osi 



2.266E-06 
1.290E-05 
O.OOOE-00 
O.OOOE-00 
5.407E-07 
1.877E-03 
4.724E-06 
1.518E-04 
O.OOOE-00 
O.OOOE-00 
2.666E-05 
5.035E-04 
6.641E-05 
7.599E-05 
O.OOOE-00 
O.OOOE-00 
4.961E-05 
O.OOOE-00 
6.550E-04 
3.445E-05 
2.349E-05 



mental composition from |Asplund et al. ( 2005)1 but the neon 
abundance from Cunha et al. (2006) and the isotopic percentage 
from Lodders| ( |2003] ). Three models were computed with X=0.7, 
Y=0.28, Z=0.02 for comparison with models of Heger & Langer 
( 2000| l. To investigate the influence of lower metallicity on the 
surface mixing an additional 12 M© model was calculated with 
X=0.736, Y=0.257, Z=0.007. 

In Table|2]some parameters of the performed models are pre- 
sented. In the first column the initial stellar mass Mi„i is listed 
and thereafter, equatorial velocity over critical velocity i/ini/^'crit, 
initial angular momentum Ji„i, MS life time th, equatorial ve- 



3. Models 

3.1. Rotation and mixing 

In Fig. [lithe diffusion coefficients in a 12 Mq model with in- 
termediate rotation velocity ((Veq) ~200 km s"') are plotted 
for three different times on the MS. When meridional circula- 
tion currents descend at the equator and ascend at the pole, i.e. 
when it transports momentum towards the centre, Dgff is drawn 
as (green) dashed line and for the opposite circulation direction 
as dotted line. 

All the models begin their evolution on the ZAMS with a 
flat angular velocity profile (Q ^constant). The profile of Q con- 
verges rapidly towards an "equilibrium" profile where the advec- 
tion of angular momentum towards the inner layers is compen- 
sated by the diffusion of angular momentum towards the outer 



layers (see Denissenkov et al. 1999 Meynet & Maeder 



20001. 



The slow expansion of the stellar outer layers during the MS 
phase, the core contraction and the effects of meridional currents 
and shear diffusion, lead to a continuous and slow change of Q. 
Very early in the evolution, a situation with two cells of 
meridional currents sets in: an inner shell which brings angular 
momentum towards the surface and an outer cell which transport 
angular momentum inwards (see Fig. fTJi). Close to the end of 
the MS even more meridional current cells appear (see Fig. fit). 
Except for a short while at the very beginning of the evolution, 
the transport of the chemical species is mainly due to Dgfi near 
the convective core and to Dshear in the outer part of the radia- 
tive zone. Dshear is thus the key parameter responsible for boron 
depletion at the surface. The nitrogen enhancement at the sur- 
face is due to the effects of both Detj and Z)sheai since nitrogen 
is enhanced in the convective core and thus must be transported 
through the whole radiative envelope. 



3.2. Evolution of surface composition 

In Fig. l2] the variations with the Lagrangian mass coordinate of 
various elements are shown. At the beginning of the evolution. 
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Fig. 2. Mass fractions of light isotopes versus mass for a 12 M© star with vini/t'crit = 0.4 at the start of hydrogen burning (a) when 
central hydrogen mass fraction is 0.48 (b) and 5 ■ 10"^ (c). In this model the mixing is strong enough to deplete boron almost 
completely at the end of the MS, lithium and beryllium are depleted at much earlier times. 



on the ZAMS, there is a thin surface layer containing Li, Be 
and B. The mass of this surface layer is between 0.3 and 2 M© 
depending on the isotope and the model considered. We see that 
the isotopes of Be and Li disappear very rapidly from the surface 
(and therefore from the whole star!). The boron isotopes are also 
rapidly depleted at the surface although on a longer timescale 
than Li. Changes of the abundances of -^He and ^"^N at the surface 
take still more time. 

The physical reasons for these different timescales associ- 
ated to the changes in surface abundances are the different tem- 
peratures required to change the abundances of these elements 
by nuclear reactions: the Li isotopes are destroyed at about 
3 ■ 10^ K, Be and B isotopes start to be depleted as soon as the 
temperature reaches about 5 and 6-10^ K respectively, ^He and 
'"^N still need higher temperatures to be destroyed/synthesized 
(in case of nitrogen) of the order of L4 ■ 10^ (for ^He) and of 
L7 ■ 10^ K (for nitrogen). In the absence of any mixing in the ra- 
diative zones, as would be the case in standard non-rotating mod- 
els, the surface abundances of these elements would not change. 
When some mixing processes are at work, as those induced by 
rotation, one expects changes of the surface abundances, more 
rapid for those elements whose abundances are changed at low 
temperatures, near the surface (such as Li, Be, and B), than for 
those which are depleted/synthesized in deeper layers (such as 
^He and '"^N). For instance, boron depletion is obtained by trans- 
port processes in a much smaller portion of the star than the one 
required to obtain nitrogen enhancement. Nitrogen indeed needs 
to be transported through the whole radiative envelope, while 
boron only needs to be transported through a small part of it. 
Therefore present models predict the existence of boron depleted 
stars with no nitrogen enrichments. From Fig.l2] we see also that 
the ratio of "B/'"B increases when evolution goes on, because 
'°B is destroyed closer to the surface than "B but this features is 
probably not observable since it occurs when both isotopes are 
already strongly depleted. 

Figure [3] shows how boron depletion occurs at the surface 
when different initial rotation velocities (and therefore angular 
momentum content) are considered. In this diagram, evolution 
proceeds from right to left. Non-rotating models would show 
non-depleted boron surface abundances during the whole MS 
phase. Only when the star is at the red supergiant stage the model 



predicts a lowering of the surface abundance in boron (loge(B n 
drops down to 1 .59). This is due to the dilution of the boron-rich 
outer layer with deeper boron depleted layers when an exter- 
nal convective zone appears. The evolution in Fig. l3] of rotating 
models is quite different with depletion of boron already dur- 
ing the MS evolution. We can see a first phase during which 
the surface velocity decreases, while no changes of the surface 
boron abundance occur. The time spent during that phase de- 
pends on the initial rotation. As a numerical example, this first 
phase lasts about 6 Myr for the 12 M© model with Wini/Vcnt = 0.1 
and about 0.6 Myr for the 12 Mq model with L'ini/L'crit = 0.6. 
During a second phase, the surface abundance decreases. The 
decrease occurs nearly at constant surface velocity in the case of 
the Wini/wcrit - 0.6 model, indicating that the mixing timescale 
is very rapid. It occurs on a longer timescale for lower initial ro- 
tation rates. Interestingly, we see that stars with a low initial ro- 
tation (below about i/ini/ycnt = 0.2) still have observable boron 
surface abundances (loge(B)>l) at the end of the MS phase. This 
is an interesting feature. Indeed boron on the surface of stars in 
the HR gap would tell us that these stars had a small rotation 
rates during the previous phases. Another important point is that 
boron depletion is very sensitive to the metallicity. Stars with 
sub-solar metallicity are more compact and undergo enhanced 
mixing. They end up with a stronger boron depletion for the 
same evolutionary stage on the MS. Our 12 M© model with half 
solar metallicity reaches a lower boron surface abundance by 
2 dex at the end of its MS life. 

Beyond the MS all models develop a convective zone at the 
surface, in which the remaining boron is diluted but not burned, 
since the temperature at the bottom of the convective zone is 
only about 6 ■ 10^ K. For the very slow rotators (wini/wcrit ^ 0.1, 
(Weq) < 50km s"' on the MS) our models predict observable 
boron abundances (loge(B)>l) even in the red supergiant phase. 

The evolution of boron and nitrogen abundances at the sur- 
face of our stellar models is drawn in Fig. |4] The 12 M© mod- 
els with i/ini/wcrit between 0.1 and 0.7 follow a similar path (see 
Fig. [4^), with the exception that the faster rotators have larger 
changes in nitrogen and boron by the end of the MS. Thus 
we see that a change in the initial velocity mainly affects the 
timescales for the changes of the surface abundances (more rapid 



' loge(X):=log(Y(X)/Y(H))+12 with Y(X) the number abundance of 
element X 
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Fig. 3. Boron versus equatorial velocity for 12 Mq models with different rotation velocities. The stars from the chosen sample with 
known rotational velocities are shown in (a) and all stars with their vsini in (b). The horizontal line indicates the minimum boron 
abundance that can still be detected. In (a) the dashed rectangle illustrates the section shown in (b). The red hatched area shows the 
evolution region for the post MS phases of our models (assuming i=90°), whereas the light red area could still be on the MS for 
models with sub-solar Z. The stars are labelled by the number given in the first column of Tables [3] and [4] Star 4 has either i/gq =14 
or 28 km s ' (see Aerts et al. 2006|l. In (b) the curves depict models with different inclination angles. 
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Fig. 4. Boron versus log(N/C) for 12 Mq models with different 
rotation velocities (a) and of 9, 12, 15, and 20 M© models with 
L'ini/i^crit =0.4 (b). In the left hand plot, the end of the MS of 
the slow rotators (yini/fcrit =0.1-0.3) is marked by a star symbol. 
The ZAMS position is indicated by a circle. 



with higher rotation rates) but not the correlation much between 
the abundances of these two elements during the MS phase. 



We also see that similar correlations are found for stars of 
different initial masses (see Fig.fflj). We can, however, note that 
the lower the initial mass, the stronger the nitrogen surface en- 
richment at a given boron abundance. This is because a given 
boron abundance is reached after a significantly longer time in 
the 9 Mq stellar model than in the 15 M© one, thus giving more 
time for changes in nitrogen in the surface layers of the 9 M© 
model. 

The depletion of boron is also correlated with changes in 
other abundances. The case of ^He is shown in Fig. p] The nu- 
clear reactions affecting ^He occur at higher temperature than 
those affecting boron and at a lower temperature than those af- 
fecting nitrogen. Thus the changes on the surface of ^He occur 
more rapidly than those of nitrogen but less rapidly than those of 
boron. 

It is interesting to look at possible correlations between sur- 
face abundances of boron and '^He since the abundances of 
these two elements can be obtained by spectroscopy for OB-type 
stars. Helium enrichments have, for instance, been obtained by 
[Lyubimkov et al.| (12004). From a theoretical point of view, one 
expects that the changes in surface helium abundance take much 
more time than changes in nitrogen. This comes from the fact 
that the abundance gradient of helium that builds up at the bor- 
der of the convective core is quite shallow with respect to the 
gradient in the abundance of nitrogen, and the stronger the gra- 
dient, the more rapid the diffusion (see e.g. Eq. 3 in |Meynet| 



et al. 2004 1. The gradient of nitrogen is steeper than the one of 



heUum because nitrogen is very rapidly enhanced in the core as 
a result of the CN cycle, while it takes much longer timescales to 
increase the central helium abundance. One consequence is that 
the present models predict that, as long as boron is observable at 
the surface on the MS, no helium enrichment is predicted. 
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In Fig. p] the boron versus '^C/'^C ratio shows the same 
property as for boron versus nitrogen; i.e., the curve is almost 
independent of the parameters velocity and stellar mass in the 
investigated parameter range. The initial value for the '^C/'^'C 
ratio is around 89 ( |Lodders| |2003[ ). The ratio continuously de- 
creases during the MS phase. On the surface, the '^C/'^C CNO- 
equilibrium value is only reached by the fastest rotators with a 
time-averaged equatorial velocity over 200 km s ' on the MS. 

A comparison with the models of'Heger & Langer| (|2000) is 
shown in Fig. IT] For the purpose of comparison, we computed 
three 12 Mg models with similar initial angular momentum and 
composition to the models of |Heger & Langer|pOOOj l. We see 
that the present models show more boron depletion at the end 
of the MS phase than the ones of Heger & Langer ("2000). The 



cause of this difference may be the different way of implement- 
ing the effects of rotation. An important difference is the way the 
advection of the angular momentum transport due to meridional 
currents is treated, as a diffusive process in the model by Heger 
& Langer] ( |2000| l and as an advective one in the present model. 
Also the counteractive effect of //-gradient on the shear diffu- 
sion is not treated in the same way in both models (see |Heger| 
& Langer| |2000[ [Meynet & Maederj |2000 ). We can also men- 
tion here that the models with solar like composition according 
to |Asplund et al. ( |2005| l, i.e. with Z=0.014, show even stronger 
surface mixing. Lower metal content makes rotational mixing 
more efficient (Maeder & Meynet 2001 1, since the stars are 



more compact. Beyond these differences, we obtain here sim- 
ilar qualitative results to Heger & Langer (2000), namely that 
boron depletion occurs much more rapidly than surface nitrogen 
enhancements. 



4. Comparison with the observations 

In Tables [3] and p] the physical properties and the surface abun- 
dances of a selected sample of stars having boron determination 
are presented. Out of the available OB-type stars with boron de- 
tection from B III line at 2065.8 A, we chose only those with Teff 
between 1 8' 00 and 29' 000 K. The B III 2065.8 A hne strength 
has a plateau ( Venn et al. 2002 1 in this temperature range, mak- 
ing the B -determination more precise. All stars of this selec- 
tion have narrow line spectra (see e.g. Proffitt & Quigley, 2001JI, 
which was a selection criterion to avoid blending effects. The 
vsini is therefore low (< 70 km s"'), meaning that either the 
equatorial velocity is low or that the star is seen nearly pole on. 

The observations of boron in young massive stars show vari- 
ations in loge(B) from 2.9 down to unobservable quantities be- 
low 1 ( [Mendel et al.][2006t|\^nn"etaLl|2002t [Proffitt & Quigley| 
2001 1. Their positions in the log(gpoi) versus log(reff) plane is 
shown in Fig. l8] For plotting the tracks, we used the polar gravity 
because it is not affected much by rotation and is a good indicator 
of the evolutionary stage of the stars. The "observed gravities", 
deduced from spectroscopy, of course do not necessarily cor- 
respond to the polar ones. There is equality between these two 
quantities when the star is slowly rotating. For the fast rotators, 
the observed gravity can be lower than the polar one if the star 
is, for instance, seen equator on (see the nice discussion of that 
topic in 'Huang & Gies"2006). Thus some points in Fig.lSlmight 
be shifted towards higher values (downwards) if it were possible 
to deduce the polar gravity from the observations; however, this 
effect is only important for very fast rotators (wini/^'crit > 0.7). 

We also see that there are stars below the ZAMS, with grav- 
ities above 4.3. Obviously, the above inclination effect cannot 
be invoked here since this effect would still push these stars to 
greater log^. But with the current uncertainties, this difference 
is not significant (see the sizes of the error bars in Fig.lSll. From 
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Fig. 7. Boron versus time for 12 Mq models with different rota- 



tion velocities, our models (solid lines) and models from Heger 
& Langer ( 2000| l (dashed lines) with comparable initial angular 
momentum J^j^ [10^^ erg s]. Both sets of models were calculated 
with Z=0.02. 



Fig. IS] we see that the majority of the observed stars have initial 
masses between 9 and 15 Mq. We also see that most of the nitro- 
gen enriched stars are found in the upper part of the MS band. 
This is consistent with the idea that these surface enrichments 
result from an evolutionary process. 

In Fig.l9]the boron abundances with respect to hydrogen are 
plotted as a function of gravity. The colour map illustrates the 
regions where our 12 Mq models with initial solar-like composi- 
tion show different values of log(N/C). The upper area represents 
log(N/C)< -0.4 (red/Ught grey) and the lower one log(N/C)> 
-0.4 (green/dark grey). We see that most of the stars are ac- 
counted for well by the models with solar composition and slow- 
to-intermediate rotation rates ((i^eq) =0-100km s"'), assuming 
that the stars with highest gravities show their initial composi- 
tion. The group of stars with loge(B) below about 1.7 are all 
more evolved stars, which is consistent with mixing processes 
occurring during the course of their evolution. In that diagram, 
stars with different N/C surface abundance ratios are plotted with 
different symbols. We see that, in agreement with models, most 
of the non-depleted boron stars show no nitrogen enrichment, 
and most of the boron depleted stars show nitrogen enrichments. 
This is indeed encouraging and can be taken as support for the 
mixing scenario. 

In Fig. [9] we see that two nitrogen-enriched stars (green dia- 
monds) appear in the red area where models predict no or small 
nitrogen enrichment. But those stars have large error bars on 
their N/C, indicating that they may be normal N-rich stars (see 
also stars in Fig.flOlwith loge(B)>2 and log(N/C)> -0.4), there- 



fore we shall not discuss them further in the present work. 

One also sees that two stars are B-depleted while showing 
no nitrogen enrichment (stars 7 and 16). Such stars can be ex- 
plained if the timescale for boron depletion is much shorter than 
the timescale for the surface nitrogen enrichment. Present mod- 



Fig. 8. Evolutionary tracks of our models in the log^pd versus 
logreff plane. For estimating the gravity of the theoretical track, 
we used the gravity at the pole. The gravities for the observed 
stars are those deduced from spectroscopy. Different kinds of 
points are used for stars with various surface nitrogen over car- 
bon ratios. The stars labelled with numbers correspond to the 
objects, discussed in the text. The typical error bars are depicted 
in the lower right hand corner. 



els, even those rotating very fast, have difficulties reproducing 
the surface abundances of the two stars belonging to this cate- 
gory. Probably here, another process than the processes studied 
in the present paper must be invoked. Binary mass transfer usu- 
ally occurs in case Bq allowing only transfer of B-depleted and 
N-enriched material. But in the less common case APTwe might 
imagine transfer of B-depleted, but not yet N-enriched material. 
In close binaries, some mixing might also occur through tidal 
mixing, a process that has not yet been explored so far whether 
from an observational point of view or from theory. Another 
possibility would be that these stars present a strong differen- 
tial rotation at the surface, triggered by some (magnetic?) brak- 
ing mechanism of the surface layers. This would in turn trig- 
ger efficient mixing through shear instabilities in the layers with 
a strong differential rotation, explaining the rapid depletion of 
boron without any significant enhancement of nitrogen. 

As just discussed, the bulk of the observations can be ex- 
plained by our tracks with an average rotation on the MS be- 
tween and 100 km s ' . Looking at Fig.|3] we see that the most 
B-depleted stars can be explained by 12 Mq stellar models with 
a velocity on the MS superior to about 90 km s"' and an inclina- 
tion inferior or equal to 30 degree^ 

What would be interesting is to know the fraction of stars 
showing a v sin / inferior to vi having a velocity v superior to a 
given limit, V2- To compute such a fraction, one has to apply the 



- Mass transfer during the transition from MS to the red-giant phase. 

^ Roche lobe overflow already during the main sequence of the donor. 

"* The inclination angle fix the position of the boron depletion curve 
on the MS at a given vsini, while the rotation velocity governs the 
amplitude of the depletion. 
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equation 



P{v sin / < vi with v > V2) = 



/ J f(v) sin advda 

V2 

J J fiv) sin advda 





, (5) 



where f{v) is the velocity distribution function, Umax the up- 
per limit of the rotation velocities, and tp the inclination angle 
between the rotation axis and the line of sight (equal to n/2 
when the axis of rotation is perpendicular to the line of sight). 
(fi = 7t/2 when v < vi and (f = arcsin(v\lv) when v > v\, 
so that V sin i < v\. We suppose an isotropic distribution of in- 
clination angles, so that the probability of having an inclination 
angle between a and a + da is proportional to sin ada. The de- 
nominator of Eq. Is] is proportional to the number of stars with 
V sin / < v\, whereas the numerator counts the subset of these 
stars with v > V2- 

For vi we adopted 50 km s ', which is appropriate for the 
investigated sample, a value equal to 100 km s"' was considered 
for V2 (see above), since this is the velocity required to repro- 
duce the strongest B-depletion. The value of i/max has been cho- 
sen as equal to 400 km s~ ' p] We assum ed a Gaussian velocity 
distribution as proposed by'Duft on et al.| ( [2006| l with the updated 
parameters, uq =225 km s"' and Au =145 km s"', from Hunter 
[eFal . (2009b). 

We find that 36% of stars with v sin / < 50 km s"' have ve- 
locities V superior to 100 km s '. Considering that our model 
12 M0 with Uini/Uciit = 0.2 only reaches the upper limits for B- 
depletion at the very end of its MS, V2 =100 km s"' might seem 
too optimistic. But with stricter limits, such as V2 =125 or even 



150 km s ', the probability is still 29% and 23%, respectively. 
This fits the 29% (9 of 31 stars) B-depleted stars in the sample. 

If we take only the sample with vsini <20 km s ' 
(vi =20 km s"'), then we get P = 0.21, i.e. 21% whereas in 
the observational sample 25% (3 out of 12 stars) are B-depleted 
(compare Fig. [3b), which is also in good agreement for the size 
of this sub-sample. We can therefore conclude that the inclina- 
tion effect is likely to play a role in the observational sample here 
and that the statistical properties fit our models of single rotating 
stars. 



In Fig. 10 we have plotted the observations of the B/H ver- 
sus N/C ratio, as well as the evolutionary tracks obtained in 
the present work. This diagram tests whether the concomitant 
changes of boron and nitrogen are reproduced by the rotating 
stellar models. In Fig. 10 1 we plotted models with uini/f ciit = 0.4 
to show the typical B-N/C relation. We see that these models 
could explain a large fraction of the observations, but they cor- 
respond to an average velocity on the MS of around 180 km s"', 
while the observational sample mainly consists of slow rotators 



{vsini <50 km s '). On the other hand, in panel b of Fig. [10 
models with lower velocities are shown. We see that they wouk 
provide a good fit to the points with boron abundances between 
0.9 and 1.5 and N/C ratios of the order of -0.2. These models 
have difficulty explaining stars that are more depleted in boron 
and more enriched in nitrogen. 



' The results only marginally depend on the upper integration limit 
as long Umax >400km s"'. 



The two stars 7 and 16 (see Fig. 10 1 present no or little sur- 
face nitrogen enrichment and are strongly boron-depleted. The 
range of values spanned by the different initial mass and initial 
velocity models barely explain these values, therefore these stars 
challenge our models even though the difference in N/C is within 
2-0" of the models at the upper limit. These two stars were also 
found to be a problem for single-star models by other authors 
dMorel et al.| [20081 IMendereFaLl [20061 [Venn et al.| [2002) 1. 
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Apart from these few difficult cases, the above comparisons 
are very encouraging on the whole and support the predictions 
of the rotating models. At this stage it is also interesting to note 
that star 7, which cannot be reproduced by the present models, 



is a spectroscopic binary with a period of 9.5 days (Pourbaix 
|et al.j |2004| l, so it might well be that this star was slowed down 
by spin-orbit coupling in the course of its MS evolution. Such 
a braking mechanism would also slow down or even stop the 
N-enrichment process. 

Are there any other indications that the changes in surface 
abundances we observe here are linked to rotation? To answer 
that question, measurements of surface velocities would be wel- 
come. At the moment only a few stars (5) have estimates of their 
surface velocity (see Table [3]). These stars are plotted in Fig. [3] 
Let us make a few comments on each of them. 

- Stars 3 and 4: for those stars, the present rotating models 
show too high surface velocity for their boron surface abun- 
dances. In other words, in Fig.pt, the theoretical MS tracks 
that go through the observed B/H ratio during the MS phase 
are too much to the right. This may indicate that mixing in 
those stars is more efficient than predicted by the present 
models^ This is supported by the high N surface content and 
that they are still on the MS. Models computed with a lower 
initial metal content may help to improve the situation (see 
yellow track and light red shaded area in Fig.l3t). Indeed, as 
can be seen looking at the track for the 12 M©, v/uait = 0.4, 
Z=0.5 Zq, the MS band extends well into the post-MS re- 
gions defined by the more metal-rich models. In that case 
there would be no difficulty at least reproducing star 3. But 
the iron content of these stars is not significantly different 
from the solar value. In general, most stars in the sample 



have Fe surface abundance consistent with the solar valuqj 
Star 4 is the most extreme case in showing B-depletion while 
rotating slowly and is therefore positioned in the shaded area 
of Fig. l3t. It cannot be explained even with a low-metallicity 
track. 
- Stars 6, 24 and 30: for these stars, the present theoretical 
models present characteristics that agree with the observa- 
tions. Looking at Fig. l3b, star 6 might be interpreted as a 
slowly rotating star at the beginning of its evolution, but its 
gravity indicates it is an evolved star that is rather at the end 
of the MS phase, if not beyond. Thus it seems that here the 
mixing was less efficient than shown by the present models. 
But the cr-error of boron is comparable to the boron deple- 
tion expected at this low rotation rate, so at the moment, this 
observation is not very constraining and is compatible with 
the present models. 

Except for star 4, present models explain surface abun- 
dances, velocities and gravities simultaneously. The number of 
stars is, however, small, and it would be interesting to obtain 
more velocity measurements to set such a conclusion on firmer 
ground. 

If we look at the projected velocities in Fig. [3t, it appears 
that all B-depleted stars are in the shaded regions, i.e. beyond 
the main sequence phase, which does not appear to be consis- 
tent with the log g values (see Fig. l9]l. This difficulty might be 
resolved by invoking the following reasons. 

1 . Initially rapidly rotating stars were suddenly slowed down at 
a given point by, for example, binary interaction. 

2. The efficiency of mixing is underestimated, so that stars 
can be both boron-depleted and rotating at low velocities. 
Rotation is measured only at the surface, while the shear 



* In case those stars were post MS stars, which is not confirmed by 
their present observed gravities (see Fig. [9j, then there would be no 
difficuUy for the present models to explain those stars. 



' Stars 5 and 8 have a loge(Fe) differing by more than 2 cr but do not 
show B-depletion. 
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turbulence, which governs the changes in the surface abun- 
dances in rotating models depends on the gradient of the ve- 
locity inside the star. Effects like magnetic breaking could 
trigger mixing in the outer layers and slow the surface down. 
3. The stars have low inclination angles (see the effect illus- 
trated in Fig[3]l. 

The last case was discussed at the start of this section and, from a 
statistical point of view, may indeed explain part of the discrep- 
ancy. Point 2 seems to be required at least for some stars (like 
7 and 16) only on the basis of their position in the boron versus 
nitrogen plane. The impact of binarity still remains difficult to 
assess in part because this effect produces similar changes in the 
surface abundances to those expected from rotation, at least for 
those stars showing boron depletion and nitrogen enrichmenr] 

A way to disentangle rotation and binary effect would be 
to look at the presence or not of a correlation between surface 
velocity and surface enrichments. In case binarity effects domi- 
nate, no correlation is expected between stellar rotation rate and 
boron depletion (Langer et al. 2010|), while rotational mixing 



culation as an advection process. This allowed us to predict the 
surface evolution of light elements in massive stars. We obtained 
the following results. 



would lead to stronger depletion for faster rotation rates (for the 
same mass, metallicity, and age). In Fig. l3] no correlation can 
be seen. This might look like favouring binary mixing. But for 
such a conclusion to be valid, one needs to correct for the incli- 
nation effects and one should be sure that stars of similar initial 
masses and ages form the bulk of the sampla^ which is not the 
case here, looking at Fig. [8] From the present data, therefore, 
we have difficulty concluding about the nature of the process at 
work. Rotation does appear to provide a reasonable process for 
explaining at least part of the stars, but binarity can also play a 
role. 

From an observational point of view, we may conclude that, 
to take a step further, the following informations would be 
needed: 

- The boron abundances on the surface of stars with a high 
V sin /. This is probably a real challenge to observe, but since 
u sin / is a lower limit to v, any high values imply a fast ro- 
tator. Unless the observed stars are at the very beginning of 
the main sequence, rotational models would predict that all 
these stars should be boron-depleted. 

- The abundances of ^He and the values of the '^C/'^C ratio 
in order to test the predictions of the present models shown 
in Figs. I5]andl6] Of course it would be interesting to have 
predictions of models invoking binary mass transfer and/or 
tidal mixing for comparison. 

- Observations in star clusters or of eclipsing binaries, since 
they would allow stars of the same age to be compared, with 
those for which some (more) precise indications on the mass 
can be obtained. 

- Asteroseismological data could probe the way rotation varies 
with depth, but it is only feasable for very slow rotators. 
Asteroseismology is also able, together with models, to pro- 
vide some information on the previous rotational history of 
the star ( Aerts 2008 1. This might offer very interesting hints 
to the processes at work here. 

5. Conclusions 

We implemented the Basel reaction network into the Geneva 
stellar evolution code (GENEC), which treats the meridional cir- 



- The boron depletion is stronger than in models of Heger 



& Langer (2000). This is probably because we account for 



the effect of meridional currents on the transport of angu- 
lar momentum as an advective process and not as a diffusive 
one, along with accounting for the effects of mean molecu- 
lar weight gradients on the efficiency of rotational diffusion, 
which is different in both models. 

We present expected correlations based on rotating stellar 
models between boron surface abundances and the surface 
abundances of ■'He, and the surface number ratios '"C/'-'C. 
Further observations will be able to check whether these cor- 
relations are present in real stars. 



^ Models with intermediate-to-fast rotation (v/v^^it ^ 0.3) show com- 
parable boron depletion (> 3 dex) to binary mass transfer. 

' Mixing efficiency depends not only on the rotation velocity alone 
but also on the age, mass, and metallicity. 



- We confirm the general conclusion obtained by'Mendel et al. 
( 2006 ) and [ Venn et aL] ( |2002 ) that rotational mixing can ac- 
count for most of the observations of the boron and nitrogen 
surface abundances. 

- We confirm the existence of challenging cases that do not fit 
well in the present scenario: stars 4 and 7, which do appear 
to present a much more rapid B-depletion than presently al- 
lowed by the models. 

- Even though our models can reproduce most of the obser- 
vations, the current uncertainties do not allow us to draw a 
firm conclusion about the questions whether our models mix 
enough and whether another surface mixing mechanism has 
to be invoked. We proposed some possibilities for further ob- 
servations that would help clarifying the situation. 
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Table 2. Model parameters 




















M,„, 


^cril 


^ini 


m 


I'cq 


P 


3He 


"He 


10b 


"B 


l^C 


»C 


14n 


[Mo] 




[1052 ergs] 


[10' yr] 


[kms-'] 








IVIass fraction 








9 


0.0 


0.00 


2.538 


(0)" 


0.720 


6.51E-05 


0.266 


7.79E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










0* 


0.720 


6.51E-05 


0.266 


7.79E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










0^ 


0.703 


4.40E-05 


0.283 


3.84E-11 


2.66E-10 


1.36E-03 


7.01E-05 


2.18E-03 


9 


0.2 


0.35 


2.632 


88(115)" 


0.720 


6.00E-05 


0.266 


5.85E-11 


5.25E-10 


2.13E-03 


5.97E-05 


8.09E-04 










76* 


0.720 


4.96E-05 


0.266 


1.14E-11 


1.38E-10 


1.76E-03 


1.03E-04 


1.20E-03 










1^ 


0.696 


3.07E-05 


0.290 


1.05E-12 


1.53E-11 


1.05E-03 


l.llE-04 


2.59E-03 


9 


0.4 


0.67 


2,663 


182(231)" 


0.720 


3.37E-05 


0.266 


2.83E-13 


7.42E-12 


1.31E-03 


1.55E-04 


1.67E-03 










163* 


0.719 


2.63E-05 


0.267 


1.04E-14 


4.68E-13 


9.48E-04 


1.56E-04 


2.15E-03 










1^ 


0.692 


1.97E-05 


0.294 


9.25E-16 


5.09E-14 


6.71E-04 


1.28E-04 


3.03E-03 


9 


0.6 


0.95 


2.692 


275 (348)" 


0.719 


2.22E-05 


0.267 


1.50E-15 


l.OlE-13 


8.25E-04 


1.53E-04 


2.28E-03 










256* 


0.717 


1.84E-05 


0.269 


1.07E-17 


1.57E-15 


5.49E-04 


1.28E-04 


2.74E-03 










2': 


0.687 


1.47E-05 


0.298 


6.87E-19 


1.48E-16 


4.04E-04 


l.OlE-04 


3.43E-03 


12 


0.0 


0.00 


1.605 


(0)" 


0.720 


6.51E-05 


0.266 


7.79E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










0* 


0.720 


6.51E-05 


0.266 


7.79E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










0^ 


0.684 


3.77E-05 


0.302 


2.30E-11 


2.35E-10 


1.33E-03 


7.15E-05 


2.42E-03 


12 


0.1 


0.31 


1.619 


44 (62)" 


0.720 


6.51E-05 


0.266 


3.05E-10 


1.89E-09 


2.28E-03 


2.92E-05 


6.60E-04 










36* 


0.720 


6.09E-05 


0.266 


9.37E-11 


7.25E-10 


2.19E-03 


5.21E-05 


7.43E-04 










0^ 


0.680 


3.25E-05 


0.306 


7.51E-12 


7.06E-11 


1.23E-03 


9.24E-05 


2.57E-03 


12 


0.2 


0.62 


1.638 


94 (123)" 


0.720 


5.33E-05 


0.266 


2.37E-11 


2.56E-10 


2.02E-03 


8.10E-05 


9.18E-04 










79* 


0.720 


4.02E-05 


0.266 


2.94E-12 


4.49E-11 


1.61E-03 


1.26E-04 


1.36E-03 










1^ 


0.678 


2.49E-05 


0.308 


2.57E-13 


4.71E-12 


1.02E-03 


1.17E-04 


2.81E-03 


12 


0.3 


0.91 


1.663 


144(184)" 


0.720 


3.63E-05 


0.266 


9.93E-13 


1.98E-11 


1.56E-03 


1.47E-04 


1.40E-03 










125* 


0.719 


2.67E-05 


0.267 


3.92E-14 


1.31E-12 


1.19E-03 


1.68E-04 


1.85E-03 










1^ 


0.677 


1.88E-05 


0.309 


3.32E-15 


1.34E-13 


8.39E-04 


1.38E-04 


2.98E-03 


12 


0.4 


1.20 


1.669 


194 (246)" 


0.719 


2.55E-05 


0.266 


2.50E-14 


1.03E-12 


1.20E-03 


1.74E-04 


1.80E-03 










173* 


0.718 


1.95E-05 


0.268 


2.97E-16 


2.43E-14 


9.05E-04 


1.73E-04 


2.23E-03 










1^ 


0.676 


1.47E-05 


0.310 


2.38E-17 


2.38E-15 


6.67E-04 


1.39E-04 


3.19E-03 


12 


0.5 


1.46 


1.676 


245 (309)" 


0.719 


1.96E-05 


0.267 


7.29E-16 


5.69E-14 


9.57E-04 


1.76E-04 


2.12E-03 










224* 


0.716 


1.56E-05 


0.269 


2.42E-18 


4.73E-16 


6.98E-04 


1.60E-04 


2.55E-03 










1^ 


0.676 


1.24E-05 


0.310 


8.87E-20 


4.36E-17 


5.33E-04 


1.29E-04 


3.35E-03 


12 


0.6 


1.68 


1.681 


295 (371)" 


0.719 


1.62E-05 


0.267 


3.51E-17 


4.47E-15 


7.85E-04 


1.68E-04 


2.36E-03 










276* 


0.715 


1.34E-05 


0.271 


2.72E-22 


1.39E-17 


5.51E-04 


1.42E-04 


2.80E-03 










1^ 


0.692 


1.18E-05 


0.293 


4.15E-24 


1.21E-18 


4.65E-04 


1.23E-04 


3.24E-03 


12 


0.7 


1.87 


1.698 


344(431)" 


0.718 


1.41E-05 


0.267 


3.19E-18 


5.71E-16 


6.81E-04 


1.63E-04 


2.50E-03 










337* 


0.713 


1.20E-05 


0.273 


7.73E-27 


5.07E-20 


4.54E-04 


1.27E-04 


2.98E-03 










2^ 


0.660 


9.22E-06 


0.326 


8.21E-29 


5.38E-22 


3.39E-04 


9.72E-05 


3.79E-03 


12 


0.4 


1.20 


1.667 


195 (248)''''' 


0.735 


1.95E-05 


0.258 


3.48E-16 


3.11E-14 


4.17E-04 


7.33E-05 


1.16E-03 










173*-'' 


0.733 


1.52E-05 


0.260 


7.39E-19 


2.67E-16 


2.83E-04 


6.38E-05 


1.38E-03 










Icd 


0.727 


1.41E-05 


0.266 


1.32E-20 


2.91E-17 


2.47E-04 


5.89E-05 


1.50E-03 


15 


0.0 


0.00 


1.128 


(0)" 


0.720 


6.51E-05 


0.266 


7.79E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










0* 


0.720 


6.51E-05 


0.266 


7.78E-10 


3.46E-09 


2.28E-03 


2.77E-05 


6.59E-04 










O' 


0.680 


3.14E-05 


0.305 


2.99E-16 


6.92E-11 


1.17E-03 


6.73E-05 


2.70E-03 


15 


0.2 


0.96 


1.197 


97 (129)" 


0.720 


4.88E-05 


0.266 


8.58E-12 


1.19E-10 


1.95E-03 


9.93E-05 


9.84E-04 










80* 


0.719 


3.56E-05 


0.267 


6.82E-13 


1.48E-11 


1.56E-03 


1.45E-04 


1.43E-03 










(T 


0.661 


2.01E-05 


0.325 


4.68E-14 


1.24E-12 


9.18E-04 


1.19E-04 


3.17E-03 


15 


0.4 


1.85 


1.209 


203 (257)" 


0.719 


2.23E-05 


0.267 


3.59E-15 


2.15E-13 


1.22E-03 


1.96E-04 


1.80E-03 










179* 


0.716 


1.69E-05 


0.270 


2.07E-17 


2.68E-15 


9.54E-04 


1.94E-04 


2.22E-03 










(f 


0.650 


1.21E-05 


0.336 


1.26E-18 


2.13E-16 


6.78E-04 


1.47E-04 


3.42E-03 


15 


0.6 


2.61 


1.224 


308 (387)" 


0.718 


1.42E-05 


0.268 


1.72E-18 


3.86E-16 


8.55E-04 


1.93E-04 


2.29E-03 










286* 


0.712 


1.14E-05 


0.274 


6.07E-28 


4.13E-20 


6.24E-04 


1.65E-04 


2.76E-03 










f 


0.642 


8.30E-06 


0.344 


4.15E-30 


2.82E-22 


4.49E-04 


1.22E-04 


3.79E-03 



" The surface velocity and mass fractions wlien tlie central liydrogen abundance X^ ~ 0.36 is given on the first line for each model. The velocity 

values in brackets are the values at ZAMS. 

* The surface velocity and mass fractions at the end of hydrogen burning (X^ ~ 10"^) is given on the second line for each model. The velocity was 

not taken exactly when the hydrogen was depleted but when it reached its minimum, i.e., a little bit earlier. 

'^ The third line for the 9 and 12 Mg models corresponds to the surface values during the RGB phase. 

"^ model with lower initial metallicity Z= 0.5-Zo. 

"^ The third line for the 15 Mg models corresponds to the surface values at the end of He-burning. 
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Table 3. Stellar parameters 



No 


Star 


Cluster 


Binary 


logs" 


'-pb 


V sin i 
[kms'] 


Ref. 


[kms-'] 


Ref. Mass'^ 

[Mo] 


1 


BD +56°576 


^Per 


Eclipse 


3.40 


22500 


17 


1 




13.5 


2 


HD886 






3.75 


22500 


10±1 


2 




10.7 


3 


HD 3360 


CasTauOBl 




3.70 


22000 


19±1 


3 


55±28 


7 10.6 


4 


HD 16582 


CasTauOBl 




3.80 


23000 


14±1 


2 


14/28 


8 10.7 


5 


HD 22951 


Per OB 


Visual 


4.40 


27900 


23 


4 




11 


6 


HD 29248 


OriOBl 




3.75 


23500 


36±3 


2 


6 


9 11.5 


7 


HD 30836 


OriOBl 


Spect. 


3.35 


21500 


43±3 


5 




13.0 


8 


HD 34816 


Ori runaway 




4.20 


28900 


35 


4 




13.5 


9 


HD 35039 


OriOBl a 


Spect. 


3.74 


20600 


4 


6 




9.0 


10 


HD 35299 


Ori OB la 




4.25 


24000 


8 


6 




9.1 


11 


HD 35337 


OriOBlc 




4.20 


23600 


15 


4 




9.0 


12 


HD 35468 


OriOBl 




3.50 


22000 


51±4 


5 




11.8 


13 


HD 36285 


OriOBlc 




4.40 


21900 


10 


4 




<9 


14 


HD 36351 


Ori OB la 


Visual 


4.16 


22000 


42 


6 




<9 


15 


HD 36430 


OriOBlc 




4.36 


19600 


26 


6 




<9 


16 


HD 36591 


Ori OB lb 


Visual 


4.00 


27000 


16±2 


5 




12.8 


17 


HD 36629 


OriOBlc 




4.35 


22300 


4 


6 




<9 


18 


HD 36959 


OriOBlc 




4.41 


24900 


11 


6 




9 


19 


HD 36960 


OriOBlc 


Visual 


4.33 


28900 


33 


4 




12 


20 


HD 37209 


OriOBlc 




4.13 


24000 


50 


6 




9.7 


21 


HD 37356 


OriOBlc 




4.13 


22400 


23 


6 




<9 


22 


HD 37481 


OriOBlc 




4.17 


23300 


67 


6 




<9 


23 


HD 37744 


Ori OB lb 




4.40 


24500 


39 


6 




9 


24 


HD 44743 






3.50 


24000 


23±2 


2 


31±5 


10 14.5 


25 


HD 46328 


Coll. 121 




3.75 


27500 


10±2 


2 




15.5 


26 


HD 50707 


Coll. 121 




3.60 


26000 


45±3 


5 




15.6 


27 


HD 52089 






3.30 


23000 


28±2 


5 




15.3 


28 


HD 111123 


Sco Cen 




3.65 


27500 


48±3 


5 




16.6 


29 


HD 205021 


CepOBl 


Spect. 


3.70 


26000 


29±2 


2 




14.5 


30 


HD 214993 


Lac OBI 




3.65 


24500 


42±4 


2 


45 


11 13.5 


31 


HD 216916 


Lac OBI 


Eclipse 


3.90 


23200 


13 


1 




9.7 



The typical 1-cr error for \ogg is between 0.15 and 0.2 (see e.g 'Morel et al.l [2008) . 
The typical 1-cr error for T^ff is about 1000 K (see e.g Morel et al. , 200" 
'^The masses of the stars are only rough estimates from the comparison of the models and the stellar positions in the logg-logre(f diagram under 



References : [l][\^nn et al. (20021; [2] Morel et al. (20061; [3] Briquet & Morel (2007); [4] Mendel et al. (20061; [5] Morel et al. ( 2008ll; [6]|Proffitt| 
|& Quigley|f2001| >; [7] ,Neiner et al., ( ,2003; ; [8] ,Aerts et aL, ( ,2006; ; [9] .Pamyatnykh et al., ( ,2004; ; [10] ,Mazumdar et al., ( ,2006; ; [ll] |Aerts|p'996^ 
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Table 4. Surface abundances 



No Star log6(C) log6(N) log(N/C) Ref. log6(Fe) Ref. log6(B) Ref. 



1 


BD +56°576 


7.84±0.18 


7.62+0.22 


-0.22+0.28 


r 


7.34+0.17 


1 


2.25+0.19 


1 


2 


HD886 


8.20±0.05 


7.58+0.11 


-0.62+0.12 


2 


7.25+0.16 


2 


2.23+0.20 


6 


3 


HD 3360 


8.16±0.08 


7.97+0.13 


-0.19+0.15 


3 


7.31+0.16 


3 


<1.21 


6 


4 


HD 16582 


8.09+0.08 


8.05+0.11 


-0.04+0.14 


2 


7.32+0.18 


2 


1.16+0.15 


1 


5 


HD 22951 


8.11±0.21 


7.69+0.32 


-0.42+0.38 


1" 


7.03+0.10 


4 


2.01+0.12 


4 


6 


HD 29248 


8.24±0.12 


7.87+0.09 


-0.37+0.15 


2 


7.36+0.19 


2 


2.45+0.40 


6 


7 


HD 30836 


8.19±0.09 


7.54+0.15 


-0.65+0.18 


5 


7.11+0.18 


5 


<1.00 


4 


8 


HD 34816 


8.17±0.07 


7.59+0.13 


-0.58+0.15 


1" 


7.12+0.11 


4 


2.17+0.13 


4 


9 


HD 35039 


8.36+0.03 


7.65+0.05 


-0.71+0.06 


1" 


7.24 


6 


2.60+0.30 


6 


10 


HD 35299 


8.19+0.19 


7.70+0.08 


-0.49+0.21 


1" 


7.19 


6 


2.89+0.30 


6 


11 


HD 35337 


8.31+0.09 


7.64+0.10 


-0.67+0.13 


r 


7.38+0.10 


4 


2.11+0.12 


4 


12 


HD 35468 


8.11+0.09 


7.90+0.16 


-0.20+0.19 


5 


7.23+0.14 


5 


<1.01 


6 


13 


HD 36285 


8.48+0.04 


7.77+0.04 


-0.71+0.06 


1" 


7.23+0.09 


4 


1.82+0.15 


4 


14 


HD 36351 


8.28+0.04 


7.68+0.03 


-0.60+0.05 


1" 


7.28 


6 


2.64+0.50 


6 


15 


HD 36430 


8.38+0.03 


7.67+0.09 


-0.71+0.09 


r 


7.54 


6 


2.53+0.50 


6 


16 


HD 36591 


8.19+0.05 


7.66+0.12 


-0.53+0.13 


5 


7.32+0.19 


5 


<1.32 


1 


17 


HD 36629 


8.32+0.03 


7.61+0.03 


-0.71+0.04 


1" 


7.39 


6 


2.47+0.30 


6 


18 


HD 36959 


8.37+0.02 


7.73+0.05 


-0.64+0.05 


1" 


7.29 


6 


2.51+0.20 


6 


19 


HD 36960 


8.39+0.01 


7.54+0.04 


-0.85+0.04 


r 


7.22+0.10 


4 


1.81+0.16 


4 


20 


HD 37209 


8.29+0.04 


7.55+0.05 


-0.74+0.06 


r 


7.32 


6 


2.48+0.20 


6 


21 


HD 37356 


8.41+0.03 


7.70+0.04 


-0.71+0.05 


1" 


7.32 


6 


2.46+0.40 


6 


22 


HD 37481 


8.39 


7.55+0.02 


-0.84 


1" 


7.40 


6 


<2.47 


6 


23 


HD 37744 


8.37+0.02 


7.80+0.04 


-0.57+0.04 


1" 


7.36 


6 


2.51+0.60 


6 


24 


HD 44743 


8.16+0.11 


7.59+0.14 


-0.57+0.18 


2 


7.17+0.19 


2 


2.76+0.20 


6 


25 


HD 46328 


8.18+0.12 


8.00+0.17 


-0.18+0.21 


2 


7.30+0.22 


2 


<1.36 


6 


26 


HD 50707 


8.18+0.10 


8.03+0.15 


-0.15+0.19 


5 


7.23+0.19 


5 


<1.50 


1 


27 


HD 52089 


8.09+0.12 


7.93+0.14 


-0.16+0.19 


5 


7.16+0.15 


5 


<1.62 


6 


28 


HD 111123 


8.04+0.10 


7.61+0.17 


-0.43+0.20 


5 


7.23+0.24 


5 


2.48+0.20 


6 


29 


HD 205021 


8.02+0.10 


7.91+0.13 


-0.11+0.17 


2 


7.24+0.23 


2 


<0.90 


1 


30 


HD 214993 


8.22+0.12 


7.64+0.18 


-0.58+0.22 


2 


7.30+0.20 


2 


2.10+0.20 


4 


31 


HD 216916 


8.17+0.13 


7.64+0.11 


-0.53+0.17 


1" 


7.66+0.14 


1 


2.31+0.15 


1 



"N and C values are from the compilation in |Venn et al.|pOO^ where most N and C values are corrected values originally from |Gies & Lambert| 

{1992) and |Cunha & Lambert|(T994l >. 

References: the numbers correspond to the same references as in TablelS] 



